function thetaDegrees = calculate_track_degree(points)
%计算轨迹段之间的角度值
[X,Y]=latlon2Cartesian_point_set(points);
%计算前后轨迹点之间轨迹段之间的夹角，以前后三个点为一组
thetaDegrees(1,1)=0;
for i=1:size(points,1)-2
    xa=X(i);
    xb=X(i+1);
    xc=X(i+2);

    ya=Y(i);
    yb=Y(i+1);
    yc=Y(i+2);

    % 计算两个向量 AB 和 BC，三个轨迹点构成两个轨迹段
    vectorAB = [xb - xa, yb - ya];
    vectorBC = [xc - xb, yc - yb];

    % 计算两个向量的点积和模长
    dotProduct = dot(vectorAB, vectorBC);
    modulusAB = norm(vectorAB);
    modulusBC = norm(vectorBC);

    % 计算两个向量之间的夹角，使用余弦定理
    cosTheta = dotProduct / (modulusAB * modulusBC);

    % 确保夹角在 [-1, 1] 范围内，避免浮点数计算误差
    cosTheta = max(min(cosTheta, 1), -1);

    % 计算夹角的弧度值
    theta = acos(cosTheta);

    % 将弧度转换为度数
    thetaDegrees(i+1,1) = theta * (180 / pi);
end
end